Length and temperature dependent crossover of charge transport across molecular junctions 



O 

<N 

bjo: 
< 



Ya-Lin Lo (BBM), h2 Shih-Jye Sun (?^±ft), 3 'E and Ying-Jer Kao (m^S) h2A ^ 

1 Department of Physics, National Taiwan University, Taipei 106, Taiwan 
2 Center of Theoretical Sciences, National Taiwan University, Taipei 106, Taiwan 
3 Department of Applied Physics, National University of Kaohsiung, Kaohsiung 811, Taiwan 
4 Center for Quantum Science and Engineering, National Taiwan University, Taipei 106, Taiwan 

(Dated: January 14, 2013) 

We study the electronic transport in a molecular junction, in which each unit is coupled to a local phonon bath, 
using the non-equilibrium Green's function method. We observe the conductance oscillates with the molecular 
chain length and the oscillation period in odd-numbered chains depends strongly on the applied bias. This 
oscillatory behavior is smeared out at the bias voltage near the phonon energy. For the phonon-free case, we 
find a crossover from tunneling to thermally activated transport as the length of the molecule increases. In the 
presence of electron-phonon interaction, the transport is thermally driven and a crossover from the thermally 
suppressed to assisted conduction is observed. 

PACS numbers: 72.60.+g, 73.23.-b, 73.63.-b, 81.07.Nb 
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I. INTRODUCTION 

The research on the molecular electronics is a fast pro- 
gressing field due to its great application potentials. Molec- 
ular junctions with single molecules or molecular assemblies 
sandwiched between two electrodes are of particular interests 
since they are potential candidates as the building blocks for 
the next generation of electronics J- One important issue as- 
sociated with the electron transport in molecular devices is 
understanding the role phonons play in the conduction mech- 
anism. In general, the electron-phonon interaction is ex- 
pected to modify transport properties when the phonon relax- 
ation time is comparable to the typical transport time of the 
junction. 2 Thanks to the rapid development of new experimen- 
tal techniques, it is now possible to study quantum transport 
phenomena in a single-molecule wire experimentally using 
break junctions, 3 STM, 4-6 and crossed-wire tunnel junctions. 7 
Vibronic contributions to the charge transport have been 
identified in the low-bias conductivity of molecular wires 
step-like features in current-voltage characteristics-^— and 
temperature-dependent conduction.— - — It has been argued 
that the transmission rate of the electron decays exponen- 
tially with length for short chains, which is a signature of 
tunneling transport. 18 For long chains, however, thermal ex- 
citations dominate at high temperature and a crossover to 
thermally activated diffusion with Ohmic behavior best de- 
scribes the electron transfer in this regime In molecu- 
lar junctions with fixed molecular length, a crossover from 
coherent to incoherent conduction, as temperature increases, 
has been observed; 21-23 in addition, strong experimental evi- 
dences demonstrated that there is a crossover from tunneling 
to hopping of the charge transport mechanism as the length 
of the molecule increases. ^r— It is not clear, however, what 
is the underlying mechanism of this crossover and how it is 
modified in the presence of electron-phonon interactions. 

In this paper, we study electrons interacting with local 
phonon baths in the molecular wires and we treat the on-site 
electron-phonon interaction perturbatively to the second or- 
der in the coupling strength. This may become important in 
molecules which couples to vibration modes provided by the 



local functional group and the charge transport characteristics 
can change significantly. Using the standard non-equilibrium 
Green's function (NEGF) method and a tight-binding model 
for the molecular chain, we study the conduction mechanism 
as a function of the length and the temperature. We observe 
a crossover of the electron conduction mechanism from tun- 
neling to thermally activated transport as the molecular length 
increases in the absence of the electron-phonon interaction. 
In the presence of the electron-phonon coupling, a crossover 
from thermally suppressed to assisted conduction is also ob- 
served. 

The paper is organized as follows. In Sec. HH we describe 
the microscopic model of the molecular wire, and the theoreti- 
cal framework to compute the transport properties in the pres- 
ence of the electron-phonon interaction. In Sec. [TTT1 we show 
our results and interpret their physics based on our proposed 
model and approximations. Finally, we conclude in Sec.HVl 



II. MODEL AND THEORETICAL METHOD 
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FIG. 1. One molecular chain, sandwiched in two leads, consists of 
N units and each unit is coupled to a local phonon bath. 

We model the molecule as a one dimensional chain of N 
units and each unit is coupled locally to a phonon bath (Fig.Q]). 
The full Hamiltonian of the molecular system with the local 
electron-phonon interaction can be written as 

H = Hi+ H c + H et + H ph + H e _ ph . (1) 

We assume the two semi-infinite one dimensional leads, and 
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the lead Hamiltonian is given by, 

CO 

Hi = e * C la C ha +VaJ2 C [a C 3^ , (2) 

aeL,R i=l 

where c\ a (c^ a ) is creation (annihilation) operator of the car- 
rier at the ith unit in the lead a = R, L , rj a is the hopping 
integral in the leads, and e a is the on-site energy in the a lead. 
The coupling between leads and the molecule is given by, 

H c = Yl V <*(°U d i + M , O) 

aEL,R,iel,N 

where the L(R) lead only couples to the first (TVth) unit of 
the molecule via the coupling constant, Vl^) , and d\ (di) is 
creation (annihilation) operator of the carrier at the ith unit of 
the middle wire. The Hamiltonian for the molecule is given 



by a single-orbital tight-binding model, 

N 

H el = £i4 d i + * d \ d 3 ' (4) 

*=1 (iJ) 

where t is the hopping integral between the nearest-neighbor 
units, 6i is the orbital energy at the ith molecular unit. 

For simplicity, we set the vibrational mode uj = ujo of the 
phonon baths to be the same and the electron-phonon interac- 
tion is given by, 

N 

H ph = y^ujpblbj , (5) 

where i is unit index and b{ (tf) is the phonon annihilation 
(creation) operator, and 

N 

H e _ ph =jJ2 d l d ^ b l + b i) , (6) 

i=l 

where 7 is the strength of the electron-phonon coupling. 

The current through an interacting region driven by a bias 
voltage is obtained by the Keldysh NEGF,— 



I= 2hJ dL ° Tr {[/i(^ML)ri - /«(w,A*fl)r*][G» - G»] + [T L - T r ]G<(lo)} 



where T ae L,R are related to the imaginary parts of the self- 
energies of leads, T aeLjR = (V a a (uj) - /2i& and 
/ a (w,/i Q ) = 1/(1 + e^ w_/ia ^) are Fermi-Dirac distributions 
in the reservoirs. 

The full (retarded, advanced and lesser) Green's functions 
are obtained by solving the Dyson equation, 

G = g + gVG , (8) 

where G is the full Green's function with the electron-phonon 
interaction, g is the non-interacting Green's function, and V 
is the interaction part. Following the notations in Ref. HH 
we solve the Dyson equation by keeping terms up to the sec- 
ond order in the electron-phonon coupling 7 to obtain the full 
Green's function, 

G = g(I + 7 G«), (9) 

where G^ contains all the information about the electron- 
phonon interaction, 

G? j = ({di(b i +bt);4))5 iij , (10) 

where we use superscript Q to indicate the presence of the 
phonon operators b and tf. ((a;b^)) represents one of the 
Keldysh Green's functions with (composite) operators a and 



(7) 
I 

b\ and their Fourier transforms are defined as 

G<H=i J dte^^a{t)), 

G&M = ~ij dte^iaitp), (11) 
for the lesser and greater Green's functions, and 

poo 

G r ab (cu) = -i / dte« u+i0+ »([a(t),bt] + ), 
Jo 

G a ab (co) = i[° dte^- M+ ^([a(t)^} + ), (12) 

J — CO 

for the retarded and advanced Green's functions. Moreover, 
in our model we consider only the on-site electron-phonon 
interaction; therefore, G Q is diagonal, and can be written as 

Gg=5 ij ((6 i + &!)+7Gg )%, (13) 

where 

G2 Q = ((di(b i + b\);(b i + b\)4)) . (14) 
The retarded Green's function, G W ' r , can be evaluated as, 
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GT "(a,) = -i / 
Jo 



dte iu>t {[d i (b i + b\y,(b i + b\)4} + ) 

(b4) 



(blk) 



uj — ujq — u> + iS uj -\- ujq — uj' -\- iS 



(hb\) 



(blk) 



UJ -\- UJq — UJ + 



iS uj — ujq — uj -\- iS 
(15) ' 



Collecting all the terms and the inverse of the full Green's 
function is given as 



numbers. 



G 1 = g 1 — T, ph , 

S pfc = 7 + b\) S id + 7 2 Gqq = H\ h 



(16) 
(17) 



where the bare Green's function g is defined as the non- 
interacting electron Green's function together with the self- 
energy corrections from the leads, 



g = {uj - H el - S 



lead ) 



(18) 



The lead self-energy, Hi ea d » describes the interaction between 
the leads and the molecular chain, and for semi-infinite one- 
dimensional leads, 30,31 



^lead — y ] 5j« j 



aeR,L 



\V n 



(19) 



<*i,n , (20) 



S< = zr a / a (o;) , 

S> = -ir a [l-/ a (a;)] , 



(21) 
(22) 



where n is 1 (TV) for a = L(i?). One remaining quantity 
required to evaluate the current is the full lesser Green's func- 
tion, G < , and it can be calculated through the G r and G a , 



G< = G r T, < G a , 

lead 1 pAi 



(23) 
(24) 



where the full lesser self-energy is broken up into contribu- 
tions from the leads and the phonon, and the former is given 
by Eq. (ETT) . For the phonon self-energy, T,^ h in Eq. ([T7l) , 
gives a shift in the energy. To evaluate the expectation value 
of the phonon operators, we perturb the phonon wavef unction 
to the first order in 7, 

!*„'.) = !"!>+£ <n f e ^ k) k), 



UJQ 



\ni) , (25) 



where |n^) and |n^) are phonon states and the energy differ- 
ence is proportional to the difference of their phonon quantum 



The second term X)^ 2 = 7 2 G gQ '<, where G gg '< is diago- 



nal with matrix elements, 



G9 Q '<(co) = i I dte^dh + &J)( &i (t) + bt(t)))(4di(t)) . 



= \b\bi) gu{u ~ wo) + (bibTj 9u(lo + w ) , 

(27) 

where (b^b) = (e^ — is the phonon occupation number. 
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FIG. 2. (Color online) Current as a function of 7 at cjo = 0.4 for 
different odd-numbered chains: N = 5 (circle), 7 (triangle), and 9 
(square) at V a = 0.1. The injection electron hits a resonance state at 
the corresponding 7 to produce the maximum. 



III. RESULTS AND DISCUSSION 

We consider the simplest case of the model, Eq. (|2]), where 
two electrodes are the identical, tjl = tjr = rj , and have 
the same coupling to the molecule, Vl = Vr = V c = 1. In 
the following, unless otherwise noticed, the energy will be in 
units of V c . The parameters used in this paper are cl,r = Q = 
0,7/ = — 5, t = —2. The symbol T indicates the temperature 
both in the reservoirs and phonon baths. The conductance is 
in units of Go = 2e 2 /h, the current is in units of Jo = 2e/ft 
, and the resistance is in units of Ro = 1/Go- The current in 
Eq. d7]) is driven by a bias voltage V a . We assume that the left 
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lead is connected to the ground, with the chemical potentials 
Hl = 0,/iR = —eV a . We obtain the Green's functions, for 
example, in Eq. (TT3T) by performing numerical integration in 
the real frequency domain, and the small positive S is set to 
0.01. 



A. Even-odd effect 

We start by considering short chains first. Figures [2] and 
[3] show the behavior of currents at small 7 for chains with 
odd and even sites at T = 0.025, and the phonon mode is at 
ujo = 0.4. This case is in the quantum limit as T/ujq <C 1. 
First we observe that odd-numbered chains shows significant 
larger current than even-numbered chains, as found in Ref . H2 
and Ref. Q3|. For even-numbered chains, in the absence of 
electron-phonon interaction (7 = 0) and at small bias, the 
chemical potential of the leads falls in the gap between the 
highest occupied molecular orbital (HOMO) and the low- 
est unoccupied molecular orbital (LUMO), thus the tunnel- 
ing current is small. On the other hand, for odd-numbered 
chains, there exists a state close to the lead chemical poten- 
tial, and resonant tunneling is possible, which produces a fi- 
nite current^2^ The electron-phonon interaction has a strong 
effect on the conductance of the odd-numbered chains. The 
molecular state near the lead chemical potential is shifted by 
the interaction, and a resonance occurs at a finite 7 (Fig. 0. 
The resonance is not observed in the even-numbered chains, 
and the conductance increases monotonically as 7 increases. 
One expects this even-odd effect should vanish when the num- 
ber of sites becomes large as the HOMO-LUMO gap de- 
creases with more atoms in the molecular wire. Moreover, this 
transport behavior depends strongly on these physical parame- 
ters, the Fermi level of the leads and the HOMO and LUMO of 
the molecules. Thus, one expects this even-odd effect would 
be tuned as these physical parameters are changed, and the 
effects are reversed. 
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FIG. 3. (Color online) Current as a function of the electron-phonon 
interaction with cjo = 0.4 for different even-numbered chains: N = 
4 (circle), 6 (triangle), and 8 (square) at V a = 0.1. 



FIG. 4. (Color online) Conductance as a function of length with 
7 = for different bias voltages V a = 0.3, 0.4, 0.5. Inset shows 
the maxima of the transmission probability are at N = 31, 61 and 
N = 41 for E = 0.4 and E = 0.3 , respectively. 



B. Length dependence 

We will focus on odd-numbered chains in the following dis- 
cussion. Figure @] shows the conductance versus the chain 
length at bias voltages V a — 0.3, 0.4, 0.5. It is clear that 
the conductance also shows oscillatory behavior as the chain 
length increases, consistent with literature results The pe- 
riod of the oscillation is roughly inverse proportional to the 
injection energy of the lead electron. This should be distin- 
guished from the oscillatory behavior in the even-odd effect as 
discussed in the previous section since this behavior can not 
be easily explained by the simple argument whether a state 
is present near the Fermi energy. The conductance maxima 
and minima can be attributed to the constructive and destruc- 
tive interferences of the outgoing and reflected wave func- 
tions, which can be better understood by solving the coupled 
Schrodinger equations 36,37 for transmission probabilities as a 
function of the electron injection energy. Details of the cal- 
culation are given in the Appendix. Inset of Fig. [4] shows the 
transmission probability as a function of the chain length with 
external energy E = 0.3, 0.4, and the resonances occur at 
roughly the same molecular length as in the Green's function 
calculation, confirming the above statement. 

Figured shows the conductance as a function of length with 
the phonon frequency ujo = 0.4, and electron-phonon cou- 
pling 7 = 0.1 at various of biases. At low bias, the electron- 
phonon interaction has little effects on the length dependence. 
As the bias increases 04=0.31), the maximum is suppressed 
and the peak position starts to shift. When the external bias 
is close to the phonon frequency (V a =036, 0.41), signature of 
the oscillation is completely washed out. Further increasing 
the bias above the phonon frequency, the oscillatory behavior 
reappears. Figure shows the total density of states for the 
molecules with different length N, 



D(E) = -Tr(ImG r (£)), 



(28) 
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FIG. 5. (Color online) Conductance as a function of length at V a =(a) 0.21,(b) 0.31, (c) 0.36, (d) 0.41, (e) 0.44, and (f) 0.51 for 7 = (circles) 
and 7 = 0.1 (triangles) with cjq = 0.4. 
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FIG. 6. (Color online) Density of states for molecular with N sites 
in (a) the phonon-free and (b) 7 = 0.1, ujq = 0.4 cases. 



which roughly corresponds to the conductance at bias V a —E 
(see Eq. (|7])). Comparing Fig. 0(a) and (b), it is clear that near 
E = coo = 0.4 the density of states is significantly modified 
by the electron-phonon interaction. Far away from this en- 



ergy, the density of states between the two cases are almost 
identical, and the oscillatory behavior reappears. Since we in 
the quantum limit where T/uj q <C 1, there are few thermally 
excited phonons present at low bias. Increasing the bias to 
larger than uo, physical phonons can be excited and inelastic 
transport process becomes available. However, these physi- 
cally excited phonons act as impurities by adding an onsite 
potential, 38 and the constructive and destructive interference 
patterns reappear. On the other hand, close to the threshold 
of the phonon frequency, one expects large fluctuations of the 
phonon number, and the phase coherence is lost and the inter- 
ference pattern disappears. 



C. Crossover of the conduction mechanism 

It has been observed in experiments that the conduction 
mechanism of conjugated molecular wires changes from tun- 
neling to hopping as the length of the molecule increases in 
oligophenyleneimine (OPI), 24 oligonaphthalenefluoreneimine 
(ONI), 39 oligo(p-phenylene ethynylene)s (OPE), 25 and 
oligophenylenetriazole (OPT)2&. At low bias voltage, tunnel- 
ing transport dominates for short chains and the conductance 
G is described by Magoga's law, 40 

G = G e _/3Ar , (29) 

where Go = 2e 2 /h, N is the molecular length, and f3 is 
the material-dependent inverse decay length. Experimentally, 
one way to distinguish the tunneling transport from the ther- 
mally activated transport is to compare the f3 values at differ- 
ent regimes. 41 For example, in the conjugated OPI molecule, 24 
a larger f3 = 3.0nm _1 is observed in the short wire, while in 
the long chain, a comparatively smaller value of /3 = 0.9nm _1 
is observed. It is argued that two different f3 values for molec- 
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ular wires of different lengths based on the same oligomers 
imply that there is a change in the conduction mechanism. 
In contrast, in the thermally activated regime, the resistance 
increases only linearly with the length and the conductance 
has a strong temperature dependence. Therefore, we use the 
following criteria to distinguish tunneling and thermally acti- 
vated transport in molecular wires:— (1) Large /3 with weak 
temperature dependence of the conductance corresponds to 
tunneling, and small (3 with strong temperature dependence 
corresponds to thermally activated transport. (2) Linear length 
dependence of the resistance corresponds to thermally ac- 
tivated transport, and exponential decay of the conductance 
corresponds to tunneling. 

Figure [7] shows the conductance versus length with fits of 
/3's in different regimes. First we discuss the case where 
the electron-phonon interaction is absent. Tunneling transport 
dominates in short chains with a larger f3 = 0.11 (per unit), 
while thermally activated transport prevails in long chains 
(N > 21) with a smaller /? = 0.04 and a linear growth of 
resistance with length (Fig.©. In addition, in the long-chain 
molecule such as N = 31 (Fig. [9] (a)), a crossover from tun- 
neling to thermally activated transport is observed as the tem- 
perature increases. At high temperatures, conductance grows 
as the temperature increases, which is a characteristic of the 
thermally activated transport conduction. The conductance is 
less sensitive to the temperature at low temperature, indicating 
the dominant transport is due to tunneling. For a short chain 
N = 5, in contrast, the conductance is temperature indepen- 
dent in the whole temperature regime, suggesting tunneling 
conduction. 




FIG. 7. (Color online) Conductance vs length at V a = 0.1. Triangles 
(squares) are for 7 = 0.1, and ujo = 0.2 (0.7). Circles are for 7 = 0. 
The slope corresponds to the f3 in Eq. (|29] >. 

Figure [lOl summarizes the conduction mechanism deter- 
mined as described above for different temperatures and chain 
lengths. It is clear that at low temperature tunneling domi- 
nates. At high temperature, a change in the transport mecha- 
nism from tunneling to thermally activated transport occurs 
as the length increases. Experiments show that the elec- 
tron transport mechanism changes abruptly at a particular 




FIG. 8. (Color online) Resistance as a function of length at V a = 
0.1V. Circles corresponds to 7 = , and triangles (squares) corre- 
sponds to 7 = 0.1 with cjo = 0.2 (0.7). Solid lines show the linear 
relation between the resistance and the chain length. 



Iength*2i"^i22 However, our calculation [Fig. [TO] (a)] shows 
that a crossover region appears between N=7 and N=19, 
where the characteristics of tunneling and thermally activated 
transport coexist. 18 Interestingly, for long chain molecules, 
there exists another transition from tunneling to thermally 
activated transport as the temperature increases. For short 
molecules, however, no such transition is observed and tun- 
neling is still the main transport mechanism at high tempera- 
ture. 

The electron-phonon interaction can significantly change 
the transport characteristics. Figure |7] shows the conduc- 
tance versus length for 7 = 0.1, with phonon frequencies 
co>o = 0.2 (triangles), 0.7(squares). The inverse decay length 
f3 decreases in the presence of the electron-phonon interac- 
tion for short chains, suggesting the coherent electron tun- 
neling becomes weaker as the interaction with phonons is 
turned on. Increasing uoq decreases the probability of the 
phonon processes, and the behavior of the conductance re- 
duces back to the non-interacting case (7 = or ujq — >> 00). 
In addition, for ujq — 0.2, the resistance is proportional to 
the length even in short chains (Fig. [8]) indicating thermally 
activated transport for all molecular sizes. Temperature de- 
pendence of the conductance [Fig. [9] (b)] shows an interest- 
ing crossover from a thermally assisted transport in the long 
molecule (N = 31) to a thermally suppressed transport in the 
short molecule (N = 5), and the transition occurs roughly at 
TV = 11. This crossover of the temperature dependent con- 
duction as length increases replaces the change of the trans- 
port mechanism in the phonon free case [Fig.[TOl(b)l. 

To better understand the temperature dependence of the 
conductance at bias V a , formally defined as 



g{v a ) = lim 

e+->0 



IJYa + e + ) - IJYa ~ e + ) 
2e+ 



(30) 



where e + is an infinitesimal positive number, we perform fur- 
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FIG. 9. (Color online) Conductance as a function of the inverse tem- 
perature for different lengths: N = 5 (circle), 11 (square), and 31 
(triangle) at V a = 0.1. ojo = 0.2 and (a) 7 = and (b) 7 = 0.1 . 



ther analysis as follows. When the electron-phonon coupling 
is weak, the first term in the current definition (Eq. (|7}) dom- 
inates and the density of state is almost the same for the two 
bias voltages separated by the infinitesimal e + , thus 



g(V a ) * lim 



duj Tr{F(K)D(^)} 



F(y a ,T,L>) = f(V a + e+ 7» - f(V a - e+ , 
i(V a ,T,u>) = f L (uj^ L )T L - f R (uj^ R )T R , 

r>(uj) = i[G r (uj)-G a (uj)} , (31) 

where TrD(cj) gives the total density of states of the molec- 
ular chain, including the effects from electron-phonon cou- 
pling; f(y a ,T,Lj) is the net electron distribution from two 
leads at temperature T and bias V a . The derivative of the elec- 
tron distribution at V a is given by lim e +^ F(V a ,T, uj)/2e + . 
We define a parameter T as 

T(K,7i,T 2 ) = J ^Tr{AF(K,T 1; T 2 ,a;)D(a;)} 



dcjAFnMAi^) + AF nn (uj)D nn (uj) 

(32) 



35 




l 


1 1 1 




1 1 1 


1 1 1 1 


| 1 1 1 1 


1 1 1 




30 




fa 1 


AAA 


A 




O 


O 


O 


= 


25 


























AAA 


A 


















AAA 


A 










: 


20 






□□□ 


□ 




O 


O 


O 










n~i □ 


□ 













- 


15 






□□□ 


□ 




O 


O 


O 










□□□ 


□ 


O 


O 


O 


O 




10 






□ □ 


□ 


O 


O 


O 


O 










EE 




O 


O 


O 


O 




5 






COO 


O 




O 


O 


O 









1 


1 1 1 






1 , , , 


1 , , , , 


1 , , 





VJT 





1 1 1 






1 1 1 1 


| 1 1 1 1 


1 1 1 






(b)~* 


A 






O 


O 






AAA 


A 






O 


O 






AAA 


A 






O 


O 






AAA 








O 


O 






AAA 






O 


O 


O 






AAA 




O 


O 


O 


O 










O 


O 


O 


O 






TTT 




O 


O 


O 


O 






TTT 




O 


O 


O 


O 






TTT 




O 


O 


O 


O 






, , 1 






1 , , , 


1 , , , , 


, , , 





VJT 

FIG. 10. (Color online) Charge transport mechanism for different 
length and inverse temperature at V a — 0.1 for (a) 7 = 0. Cir- 
cles indicate tunneling transport, triangles correspond to thermally 
activated transport, and squares represent the crossover regime with 
the mixture of tunneling and thermally activated transport character- 
istics (see text), and (b) for 7 = 0.1 and cjo = 0.2. The up and 
down triangles indicate thermally assisted and suppressed conduc- 
tion, respectively, and stars represent a transition between these two 
behaviors. 



where the difference between distributions at two differ- 
ent temperatures is AF(V a ,T 1 ,T 2 ,uj) = F(V a ,T u uj) - 
F(Va,T2,co>), and has only two non-zero matrix elements 
AFn, and AFnn due to the coupling to the left lead Tf x 
and to the right lead respectively. 

The physics of a thermally activated or suppressed trans- 
port is thus dictated by T. In the case T x > T 2 , T(V a ,T 1 ,T 2 ) 
is positive for thermally activated transport , and negative for 
thermally suppressed transport. Since the density of states has 
small temperature dependence, AF(uj) would play the major 
role on the temperature behavior. Assuming the density of 
states at sites 1 and N are approximately equal due to sym- 
metry, we can rewrite Eq. (l32t as 



T(K,Ti,T 2 )« J duj |AFn(u;) + AF NN (u) D n (cj) 
= J dujAF(uj)Du(uj), 
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where F = Fu + Fnn- Inset of Fig. [TT] shows the 
AF(Va,Ti = 300K,T 2 = 70K,cj) as a function of energy, 
We notice AF is negative near uo — —V a , and becomes posi- 
tive away from this region. Figure [TT] shows the integrand of 
Eq. (l32l) . For N = 5, there is more negative area than the pos- 
itive one, which gives T < 0, and the transport is thermally 
suppressed. On the other hand, for the N = 31 chain, there 
is more positive area than the negative one, resulting in a ther- 
mally activated transport. In summary, the temperature depen- 
dence of the conductance comes principally from AF(u), and 
the characteristics of the molecules and the electron-phonon 
effects enter through the density of states D(uu). For N = 5, 
the weak electron-phonon interaction results in a shift of the 
density of states toward u = —V a . This picks up more neg- 
ative contribution in AF, resulting in a thermally suppressed 
transport. 




0.2 



FIG. 1 1 . (Color online) The integrand in Eq. (l32l as a function of 
energy uj for N = 5 (circle) and N = 31 (triangle) at Va = 0.1, and 
Ti=300K, T2=70K. Inset: Temperature difference of the effective 
electron distributions of the leads AF(cj). See text for details. 



in long chains. For a long molecule of fixed length, there 
exists also a crossover from tunneling at low temperature to 
thermally activated transport at high temperature. In the pres- 
ence of the electron-phonon interaction, the coherent tunnel- 
ing is interrupted by the phonon and the length crossover dis- 
appears. The thermally activated transport becomes the main 
transport mechanism. At high temperature a crossover of tem- 
perature dependent conduction from the thermally suppressed 
to assisted transport as the length increases, in place of the 
crossover in the transport mechanism in the phonon free case. 

Without resorting to the ab initio machinery, we are able 
to extract essential physics in these molecular junctions us- 
ing a simple model, and the results qualitatively agree with a 
wide range of experiments on molecular wires. The approach 
in this paper can be generalized to study molecules with site- 
dependent coupling strengths and phonon energies. This may 
be important to understand the mechanical roles of the func- 
tion groups. In addition, the results can be extended to the 
spin dependent case to study the properties of the magnetic 
tunneling junctions. Further improvement can also be made 
by self-consistently calculating the charge density at each site. 
These questions will be addressed in a future study. 
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IV. CONCLUSION 

In conclusion, using non-equilibrium Green's function 
method, we study the length and temperature dependences 
of electron transport in molecular junctions. We model the 
molecular wire using a simple tight-binding model, with each 
site coupled locally to a phonon bath. We treat the electron- 
phonon interaction perturbatively up to second order in the 
electron-phonon coupling. 

We find different behaviors of conductance in short chains. 
Larger conductance is observed in odd-numbered chains, and 
a resonance is present when the electron-phonon coupling in- 
creases. On the other hand, for even-numbered chains, low 
conductance is observed with a monotonic increase as the 
electron-phonon coupling increases. In the absence of the 
electron-phonon interaction, we observe a length dependent 
crossover of the conduction mechanism from tunneling to 
thermally activated transport at high temperature. We find 
that at high temperatures, tunneling transport dominates in 
short chains, while thermally activated transport dominates 



Appendix A 

In this appendix, we provide an alternative method to study 
the transport properties of molecular chains by solving cou- 
pled Schodinger equations in the tight-binding limit for the 
on-site wave functions. For illustrative purposes, we consider 
the molecular transport without the electron-phonon interac- 
tion, and the full Hamiltonian can be written as, 

H = H l + H c + H eU (Al) 

where the Hamiltonian of two leads H\ , the molecular elec- 
tronic system H e i , and the coupling between them H c are 
described by Eq. © to ©. We solve the site -dependent 
Schrodinger equations Efa = Y^j Hij 0j w * m me tight- 
binding model to obtain the transmission and reflection am- 
plitudes. In the leads, the dispersion relation E — e a = 
2r] a cos(k a a a ) for a = L, R is employed, where e a includes 
both the on-site energy and the chemical potential in a lead, 
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a a is lattice constant, and ki (k r ) is incoming (outgoing) wave 
vector. We use plane waves as the wave functions in leads, 



o\e 
a 3 e l 



a 2 e 



(A2) 
(A3) 



where in the left lead the total wave function </>f consists of 
two components of the incoming and reflected waves, and the 
transmitted wave in the right lead. The transition and reflec- 
tion probabilities are defined as 



T = c 
R = c 



(A4) 
(A5) 



where c is a normalization constant given by |T| 2 + \R\ 2 
The current across the molecular junction is given by 



1. 



dET(E)[f L (E,iJL L )-f R (E,VL R )] 



(A6) 



After solving the Schrodinger equations, the transition am- 
plitude as/ai and the reflection amplitude 0,2/0,1 are obtained 

by, 



03 = 2 VL V L V R cos{kiai)e ik ^e- ik ^i /^A 

01 (E - e R - rj R e ik r»r)[vl - (E - e L - r] L e ik ' a ')(E - ei - tD 2 )} \ J ' 
^2 = (E-e L - VL e- ik ' a ')(E - Cl - tD 2 ) - V? 2ifciX , 
01 Vl - (E - e L - r] L e ik ' a ')(E - a - tD 2 ) 



The ratio of wave functions at the iVth and the first site is 
obtained by solving the site-dependent Schrodinger equations 
of the molecular chain, which gives 



N 



N 



where the Green's function on the nth site is given by 

t 



D n 



E-e n ' 



(A9) 



(A10) 



which propagates the wave function from the n — 1th to the 
nth site as <\> c n = D n (\) c n _ x . Moreover, the molecular chain is 
influenced by leads through the weak coupling V a , and these 
effects are included in the self-energies. We replace the on-site 



energies at the ends of the chain (n = 1 , N) with effective on- 
site energies, which include the self-energies from the leads. 
For example, the Schrodinger equation at site N is given by, 



Ecf) c N = t<i>\ 



'N-l 



N 



V R $ 



(All) 



and the Green's function Dn = t/(E — ejv)is obtained with 
the effective on-site energy cn = £n + Vr^i I ^n > wnere me 
ratio I (j) c N is obtained by solving the Schrodinger equation 
of the right lead, 



V R 



E-€r- r] R e 



(A12) 



On the other hand, the effective on-site energy cn can also 
be described by the self-energy from the right lead, cn = 
€n + ^r , and T,r can be obtained from Eq. (l20l) . 
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